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1.       INTRODUCTION 

Solid  body  equations    to   calculate   the   trajectory  of  a  projectile 
which  is   flying  with  thrust  of  ramjet  under  hypersonic  Mach  number  are 
formulated.      The  projectile  is   gun  launched.   This  projectile  has   an  axi- 
symmetric   configuration  and  is   fin  stabilized  with   four  pop-out  wings. 
The  projectile  has   canards    to  control  attitude.      Ideally   the  projectile 
has  no  spin;   however,   due   to  effects   from  miscellaneous    forces   and  moments 
which  work  on  the  body,   it  may  acquire  spin  motion. 

In  order  to  calculate  the  trajectory  considering  such  conditions, 
solid  body  equations  with  6    degrees   of  freedom  are  induced.      One  of  the 
main  problems  which  will  occur  during  flight  motion  of   the  projectile  is 
trajectory  errors  which  are   caused  by  numerous  origins  especially  due 
to  the  wind. 

Total  range  of  these  projectiles   are   rather  great.      Consequently, 
perturbing  influences   from  these  origins   are  important  to  notice  when 
long  range   trajectories   are   calculated.      Some  of   the  primary   origins   of 
error  are  mass    unbalance   due   to  impreciseness   of  construction   and  from 
shock  of  launching,   misalignment  of  jet,  misalignment  of  intake  air  flow, 
and  so  on. 

In  order   to  minimize   these  effects,    it  is   useful  to  set   the  value   of 
thrust  of   the   ramjet  equal  and  opposite   to  that  of  aerodynamic   drag. 
The    trajectory  is   equivalent    to  projectile   flight   in  a  vacuum  and  is 
known  as    thrust-equal-drag   trajectory.      The  effect   of  wind,  when  the 
direction  of  wind  is   the  same  as   that  of  the   trajectory,    can  be  readily 
calculated.      In   case   of  a  side  wind,    the  projectile  will  be  carried 
laterally. 


la  order   to  cancel   the  effect  of  a  side  wind  which  will  carry   the 
projectile  laterally   and  to  restore   the  planar  motion  of   the   trajectory, 
active   control  can  be  used.      For  example,   canards   can  be  employed.      Due 
to   active   aerodynamic  control,   total  range  will  be  affected.      Active 
control  modulates   aerodynamic  drag.      The  boundary  and  initial  conditions 
will  influence  whether  total  range  will  be  expanded  or  shortened  by   the 
effect  of  canards. 

In   this  paper,   effects   due   to  winds    are   considered;  however,   pertur- 
bations of  other  origins  are   concentrated  in  one   term  in  the  equations   of 
motion  as    constants.   In  order  to  calculate   these  data  by  electronic 
computer,    flow   diagrams   for  FORTRAN  are  indicated. 
2.     Coordinate  Systems 

In  order  to   construct  solid  body  equations,    coordinate  systems   are 

(2) 

considered  as    follows: 

(1)  Ground  axis   coordinate  system  (0-xyz) . 

The   origin  is   on  the  center  of   gravity  of   the  projectile   at   the 
time  of  departure  of  projectile.      The  x  axis   is    toward  the 
direction  of  horizontal  line  of  shot.      The  y   axis   is  vertically 
upright,   and  the  z  axis  is  horizontally  rectangular  in  a  right 
hand  sense. 

(2)  Velocity  axis   coordinate  system  (0-x  y  z   ). 

The   origin  is   at  the   center  of  gravity,    CG,   of   the  projectile. 
The   x     axis   is    along   the  velocity  vector  of   the  projectile,   i.e., 
x       is   tangent  to  the  trajectory.     The  transformation   to   create 
the   (x   ,y    ,z   )    coordinate  system  from  the    (x,y,z)    coordinate 
system  is,    first,   a  rotation  about   the  z  axis  by   an   angle   0. 
A  new  coordinate  system  is   formed  and  is   identified  as    the   (x   ,y    ,.z) 
coordinate  system. 


The  velocity  vector,   v  and  x\    are  located  in   the  x^-z  plane. 
Now  rotate  the   (x-.,y..,z)    coordinate  system  about   the  y.    axis 
until   the  x..    axis   is  parallel  with  v.    The  angle  of  rotation   about 
the  y.    axis  is  Q.      The   resulting  coordinate  system  is    the 
(x   ,   y    ,    z   )    system.      Refer   to   Figure  1   for  an   illustration  of 
the  axes. 


►  X 


Figure  1.   Illustration  of  geometry  for  the  velocity  axis  coordinate 
system. 


(3)  N on rot a ting  body  axis   coordinate  system  (0  -   x_y  z_) . 

The   origin  is   at  the  center  of   gravity  of  projectile:      the 
x^   axis    coincides  with   the  longitudinal  axis   of  synmetry.      The 
y_    and  z_    axis   are   transformed  as    the  same  manner  as    that  of 
velocity  axis   coordinate  system  by  <J>  instead  of  9,   and  by  5  in- 
stead of  a. 

(4)  Wind  axis   coordinate  system  (0  -   x     y     z     )♦     The  system  moves  with 

J  varva  va  J 

relative  wind  speed  similar  to  velocity  axis  coordinate  system. 

The  relevant  angles   are    (9    .0   )    instead  of   (9,a)  with   the   x 

°  aJ   a  *  va 

axis   along  the  relative  velocity  vector. 
3,.    Equations 

Solid  body  equations  with  six  degrees   of  freedom  are   developed  by 
the  method  of  coordinate   transformation  in  order  to  calculate   the   trajectory 
of   the  projectile     with  ramjet  propulsion  as    follows: 

x  =  v  cos  9   cos  a                                                                         (1) 

y  =  v  sin  9    cos  a                                                                               (2) 

2  =  -v  sin  a  (3) 

v  =  -  F  (4) 


m     xv 


(5) 


mv  cos   a     yv 


a  =  -  —  F  (6) 

mv     zv  v    ' 

p  =  r   /   cos   E,  (7) 

k  =  q  (8) 

n   =  r   tan  £  (9) 

b   =  MxB/Ix  (10) 

q   =    (1  -   WDpr  +  MyB/I  (11) 

r  =  -(1   -   Ix/I)pq  +  MzB/I  (12) 


F       =  -mg  sin  8   cos  a  +  T  cos    (Cp    -  8)   -   hpv2SC  (13) 

2 
F       =  -n»g  cos   8   -  Ta  +  hpv  SC     a  +  C   cos  tt  (14) 

y "  u  £ 

a 
2 
F        =  -mg  sin  8   sin  a  -  T  8+  ^pv  SC      3  +  C   sinir  (15) 

Mxfi  =   0  a  (16) 

MyB  =  ?1PV  SdCm  2  -  Trt  3  +  M£  cos  tt  (17) 

2    a 

M,«  =  ^pv  SdC   a  -  Tr  a  +  M  sin  tt  (18) 

2d  m        t      £ 

a 

a  =  tan  "1  cos  gsin(<?-6) 

cosa  cosg  cos(<^-8)  +  sina  sing 


ft  =  tan~*   sina  cosg  cos  (</M3)  -  cosa  sing 
cosa  cosg  cos(<^-9)  +  sina  sing 


(20) 


For  a  projectile  with  ramjet  propulsion,  eq.  (16)  is  true  usually.   At 
launch,  p  =  0,  and  there  is  no  initial  spin  so  that  p  =  0.   For  case  of 
p  =  0,  equations  (11)  and  (12)  become 

q  »  M  B/I  (11)' 

r  -  Mzfi/I  (12)' 

Furthermore,  if  thrust  has  a  value  near  to  drag,  then  dispersion  becomes 
a  minimum  so  that  the  last  terms  of  equations  (14),  (15),  (17)   and 
(18)  are  deleted  from  these  equations. 
Effect  of  Wind 

When  a  wind  blows,  a  wind  axis  coordinate  system  is  considered. 
Using  these  coordinates,  the  x  axis  is  along  the  relative  velocity  between 
projectile  and  air.   The  motion  of  the  projectile  can  be  calculated  the 
same  way  as  in  the  case  of  no  wind.   For  this  case,  we  use  a  subscript  (a) 
which  is  applied  to  the  notations  of  the  case  of  no  wind. 

In  order  to  obtain  velocities  which  were  observed  from  the  ground, 
we  must  add  distances  transported  by  wind  during  the  time  interval  to 
values  calculated  on  the  wind  axis  coordinate.   The  relations  between 
velocities  on  wind  axis  coordinate  and  velocity  axis  coordinate  are  as 


v  cos  9  cos  a  =  v  cos  9  cos  a  -  w        (21) 
a      a      a  x 

v  sin  9   cos  a  =  v  sin  9  cos  a  -  w         (22) 
a      a      a  y 

v  sin  Y  =  v  sin  a  +  w  (23) 

a      a  z 

From  these  aquations,  we  can  calculate  v  ,  9  ,  a  knowing  data  without 

a   a   a 

wind,  v,  9,  a  and  the  component  of  wind  relative  to  the  ground  www. 

In  the  case  of  computer  calculation  by  a  step-by-step  method  we 
can  use  (21),  (22),  (23)  to  get  initial  value  of  v,  9,  a      in  each  suc- 
cessive step.   These  values  of  v,  9,  a  are  relative  to  the  ground  axis 
coordinates,  so  that  the  trajectory  is  obtained  from  them. 
4.  Effect  of  Canard 

From  mentioned  above,  it  is  clear  that  projectile  is  not  within  the 
plane  of  launching  but  flies  in  a  different  direction  from  the  first 
schedule  when  wind  is  blowing.   In  order  to  cancel  the  miss  distance, 
we  can  use  a  canard  giving  some  suitable  trim  angle.   Trim  angle  will  be 
given  on  body  axis  coordinate  and  whose  effect  will  act  as  sinusoidal 
form  such  as  V  sin  y  where  y  is  factor  angle  measured  from  y  axis  right 
wise  around  x^  axis.   In  order  to  cancel  the  effect  of  sidewise  dispersion 

•  « 

z  =  0soa  =  0,a_=*0  must  be  satisfied, 
therefore 

F   =0  (24) 

In  order  to  demand  (24)  we  must  add  a  term  for  the  canard  on  the  right 

terms  of  eq.  (15) 

2 
-mg  sin  9  sin  a  -  TS  +  %)v  SC   3  +  C  sin  tt  =  v  sin  y   (25) 

J_i  £ 

a 
There  are  many  factors  to  give  effect  to  the  left  side  terms  of  this 

equation,  so  it  is  difficult  to  discuss  analytically    However,  the 

value  of  the  terms  are  calculated  on  each  step  so  that  we  can  decide  a 

suitable  value  of  canard  angle  for  each  step  of  the  calculation. 
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5.      Program  Flowcharts 

In  order  to  calculate   these  equations,   a  computer  program  will 

(31 
be  built.     It   is   illustrated  by  program  flowcharts.     This   computer  program 

consists  of  a  main  program  and  subroutines. 

Inputs  which  are  used  for  calculation  are   ready   and  stored  in  incre- 
ment address  by  orders  in  the  program  to  use  for  calculation.      Data 
calculated  are  provided  printed  on  sheets   and  so  on.      Program  flowcharts 
of  main  program  and  subroutines   are  shown  on   figures,   and  they   are  explained 
as   follows; 

a.  Main  Program 

Main  program  has   a  job   to  store  and  supply  data  from  incre- 
ment address  of  computer  and  to  make  initial  condition  in  shape 
suitable   for   calculation.      It   gives    these   data  of  subroutines    to   com- 
pute and  to  get  results   of  calculation.      Data  returned  from  sub- 
routine  are  printed  periodically.      Calculation  is   carried  by  one 
step  of  variables. 

These  data  are   accumulated  step-by-step   until   they   come 
to  some  objective  point.      When  the   objective  point  is  satisfied, 
the   calculation  is  stopped.      The  point  where  the  projectile 
touched   the   ground,   point  of  print  and  points    to   change   initial 
conditions   are   detected  and  ordered  to   computer.      It    contains  FORMAT 
to  print  these  results, 
b.      Subroutine  COMPUTE 

Da.  order  to  calculate  eq.    (1)   *>  (12)    there  need  data  of  eq.    (13) 
^   (20)   which  are  in  the  state  suitable   for  input.      In  this  subroutine 
these  data  are  prepared  using   three  subroutines   of  ATOM,    FORCE  and 
MOMEN.      These  data  are  sent  to  subroutine   RKG.      There  eq.    (1)   ^   (12) 


are   calculated  and  results   are  sent  back  to  this   routine.      Additional 
data  are   calculated  here.      These  data  are  sent  back   to  main  program, 
t.      Subroutine  RKG 

This  subroutine  has   a  job   to  calculate  each  of   the  dif- 
ferential equations   in  a  Runga-Kutta-Gill  method.      Suitable   data 
for  input  are  sent  from  subroutine  COMPUTE.      This   subroutine   calls 
differential  equations   one  by  one   from  subroutine   DANDO. 

d.  Subroutine  DANDO 

This  subroutine  stores   each  equations  whwch   composes    the   tra- 
jectory equations  in  state  of  FORTRAN  form. 

e .  ATOM 

This  subroutine  provides  air  temperature  and  air  density  for 
the  standard  state  of  heights  from  sea  level.  Air  temperature  is 
used   to   calculate  Mach  number. 

f.  FORCE 

Data  of  forces   are  stored  in  increment  address.      On   this   sub- 
routine  forces    that  are  arranged  suitably   to   the  state  of  input  of 
the  instance  are  selected  and  prepared  for  supply.      This   subroutine 
has    a  function  to  make  linear  interpolation. 

g.  MO  MEN 

Data  of  moments   are  stored  in  increment  address.      Moments    that 
are   arranged  suitably   to  the  state  of  input  of   the  instance   are 
selected  and  prepared   to  supply.      This   subroutine  has   a  function 
to  make   linear  interpolation. 


Program  Flow   Charts 
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(C)      Subroutine   RKG 
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(d) 
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(e)  Subroutine  ATOM 
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(f)  Subroutine  FORCE 
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(g)  Subroutine  MOMEN 
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6.  Conclusions 

Solid  body  equations  with  6-degree-of-f reedom  are   formulated  in  order 
to   calculate   the   trajectory  of  a  projectile  which  is    flying  with  the  thrust 
of  a  ramjet. 

In  order  to   calculate  by  electronic  digital  computer;  program   flow 
charts   are   developed.      FORTRAN  programs  will  be  made   for  each  subroutine  and 
program.      Trajectory  errors   due   to   the  effect  of  the  wind  are   considered 
mainly.    Dispersion  due   to  the  wind  will  be  calculated  by  these  programs. 
If  there  need  to   cancel  the  effect  of  a  side  wind,   trim  angle  of  canard  will 
be  given  step  by  step  in  the  procedure   of   calculation. 
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